Systematic Errors in Sunyaev-Zeldovich Surveys of Galaxy Cluster Velocities 
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Galaxy cluster surveys compiled via the Sunyaev-Zeldovich Effect have the potential to place 
strong constraints on cosmology, and in particular the nature of dark energy. Here we consider 
cluster velocity surveys using kinetic Sunyaev-Zeldovich measurements. Cluster velocities closely 
trace the large-scale velocity field independent of cluster mass; we demonstrate that two useful 
cluster velocity statistics are nearly independent of cluster mass, in marked contrast to cluster 
number count statistics. On the other hand, cluster velocity determinations from three-band ob- 
servations of Sunyaev-Zeldovich distortions can require additional cluster data or assumptions, and 
are complicated by microwave emission from dusty galaxies and radio sources, which may be cor- 
QQ ■ related with clusters. Systematic errors in velocity due to these factors can give substantial biases 

' in determination of dark energy parameters, although large cluster velocity surveys will contain 

, enough information that the errors can be modeled using the data itself, with little degradation in 

. cosmological constraints. 
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I. INTRODUCTION 

CIh' The formation of galaxy clusters is a sensitive tracer of structure growth in the universe. It is widely appreciated 
that the number of clusters larger than a given mass as a function of redshift has the potential to place strong 
^ . constraints on dark energy properties: cluster number counts are one of four techniques considered in detail by 
' the Dark Energy Task Force. Currently operating Sunyaev-Zeldovich experiments P, Q will soon produce large 
. catalogs of galaxy clusters to all redshifts, potentially fulfilling this promise. But the sensitivity of cluster number 
counts as a probe of cosmology also means that systematic errors in these number counts can substantially bias 
[ any cosmological constraints 0, 0, H H, Q- Including information on the spatial distribution of galaxy clustering, 
^ ■ termed "self-calibration," can help characterize and compensate for systematic errors in number counts [1,1 [13, but 
' these techniques require demanding further observations and rely on specific simplifying assumptions about cluster 
0^ . properties. 

\l ' We have recently explored an alternative use of galaxy clusters to constrain dark energy, via their peculiar velocities 
. The kinematic Sunyaev-Zeldovich Effect measures an effective Doppler shift of clusters, proportional to both 

' the cluster line-of-sight peculiar velocity and its optical depth, so measurement of this small blackbody distortion in 
[ the microwave background radiation holds the promise of directly measuring cluster velocities with respect to the 
25 ' rest frame of the microwave radiation. Cluster velocities trace the large-scale velocity field arising from structure 
formation in the universe, and their velocities are expected to be only weakly dependent on cluster mass. Therefore 
we expect that dark energy constraints based on cluster velocities will be far less sensitive to systematic errors in 
estimating the mass limit of any particular cluster catalog. We verify this expectation here. While changing a cluster 
catalog mass cutoff by 20% can change the total number of clusters by a factor of two, it only changes cluster velocity 
' statistics by a few percent. The bias on cosmological parameters from uncertainties in cluster selection will be much 
milder for cluster velocity statistics than cluster number counts, and the power of velocities to constrain dark energy 
is significant, even for modest velocity errors as large as 500 km/sec (Tsj . 

The kinematic SZ signal measures the total cluster baryon momentum, not the cluster velocities directly. Extracting 
the velocities from SZ measurements will require additional data to estimate the cluster baryon mass [14]; cluster X-ray 
temperatures are one likely route, while a suspected tight correlation between thermal SZ flux and gas temperature 
is another. But these measurements have potential systematic errors of their own; previous studies [isl . [l6| have 
shown that X-ray temperature systematically overestimates electron temperature by 20% to 40% (this particular 
systematic difference arises because while X-ray temperature is luminosity weighted, electron temperature is mass 
weighted), while correlations seen in numerical simulations may not incorporate all of the relevant physical effects 
in real galaxy clusters. Observing the SZ signal will also be complicated by point source contamination, which may 
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induce a different systematic error. Using simple models for these errors, we show that they give potentially significant 
biases to cosmological parameters if not properly accounted for. 

In this work, we focus on two particular galaxy cluster velocity statistics: the correlation function of the velocity 
components perpendicular to the line connecting a cluster pair, {viVj) ^ (r), and the mean pairwise streaming velocity 
Vij{r); each is a function of the separation between two galaxy clusters r and redshift z. Both of these statistics 
were considered as probes of dark energy in Ref. [l3l |. (We drop the perpendicular subscript from the correlation 
function for convenience) . Sections [iTl and IIIII consider the systematic errors arising from uncertainty in the cluster 
mass selection function, and from systematic errors in velocity estimates due to misestimates in galaxy cluster physics 
or contamination by foreground emission. Section IIVI then computes the resulting biases in determining dark energy 
parameters for each source of error, while Section |V] briefly considers self-calibration techniques in the context of 
cluster velocities. 



II. SYSTEMATIC ERRORS FROM MASS MISESTIMATES 

In this section, we address the issue of mass selection in the context of cluster velocities. The velocity statistics 
depend on M^i^ through the normalization term in theoretical halo models. In order to model the effect of mass 
selection, we assume that for a given survey, cluster masses are all mis-estimated by a constant fraction. This leads 
to a corresponding difference in the inferred cluster mass threshold for number statistics and the resulting systematic 
bias in dark energy parameters studied in Ref. [3]. However, clusters of any mass generally trace the large-scale 
velocity field, so biased cluster mass estimates should have little effect on cluster velocity statistics. The two statistics 
(viVj) j_ (r) and Vij(r) can both be estimated accurately with analytic approximations based on the halo model and 
on nonlinear perturbation theory [13, [3 [11] ; a summary of these approximations is given in Ref. [13| ■ 

To approximate the effect of cluster mass mis-estimates, assume that for a given sample of galaxy clusters detected 
via the SZ effect, all of the inferred masses are off by 20%. We compute the velocity statistics for clusters with both 
the actual mass cutoff and the inferred one using the halo model, and find that this large change in mass selection 
has minimal effect: a 20% offset in minimum mass estimate gives only 2% to 4% change for both velocity statistics. 
Fig. [T] displays the difference. Therefore, even if the mass determination of clusters is uncertain at this level, it will 
result in only small changes in the underlying cosmological models selected by the data. We emphasize that this is in 
marked contrast to the case for cluster number counts. Although self-calibration techniques provide a possible remedy 
to the cluster count mass-selection bias [^, [l^ , it requires at minimum a determination of the scatter and bias in 
cluster photometric redshifts to better than 0.03 and 0.003 respectively in order for self-calibration to work [20]. The 
evolution of cluster properties with redshift also must be of an assumed form. 

III. SYSTEMATIC ERRORS FROM VELOCITY MISESTIMATES 

Aside from errors in the inferred cluster mass, which the previous section shows has little effect on cluster velocity 
statistics, the cluster velocities themselves are also subject to systematic errors. The state of some small volume of gas 
is characterized by its temperature, density, and bulk velocity. Measurements of the SZ distortions of the radiation 
passing through this gas at three frequencies often have a physical degeneracy which enable measurement of only 
two of these gas quantities [3, [Ml, [13]; in particular, this is true for the ACT frequency bands at 145, 220, and 
280 GHz. To get the third (typically the gas velocity, the quantity of interest here), additional information must be 
obtained. The most convenient source is a direct determination of gas temperature, either through X-ray observations 
or through a theoretical correlation of gas temperature and total thermal SZ distortion. In addition, extracting the 
SZ signal accurately in the presence of foreground emission, particularly from infrared point sources [23l [23. |25|. can 
also lead to systematic errors in inferred velocities if the point sources are not adequately characterized. 

Typically the X-ray temperature Tx of galaxy clusters differs appreciably from the electron gas temperature Tg. 
Numerical simulations [l^ [2^ indicates that using Tx as a proxy for leads to over-estimation of peculiar velocities 
inferred from SZ measurements by 10 to 40%. This bias is because the estimate of peculiar velocity gets weighted 
by the ratio of pressure- weighted temperature to mass- weighted temperature (l6j . In order to quantify the effect of 
the difference between Tx and Tg in the estimation of peculiar velocity, we use the relation {v)e^ oc {Te)et, for an 
unresolved cluster [l^. The quantities (..) indicates the average quantities over the beam size 9b\ for brevity we drop 
the average symbols and assume beam-averaged quantities in this Section. We assume the X-ray temperature Tx 
and the electron temperature T^ are simply related through a linear relation Te = a -t- hTx- Numerical studies find 
(a, h) ~ (0.17 keV, 0.69) when averaged within the virial radius for a cluster and (0.18 keV, 0.53) when averaged within 
three times the virial radius flS*] . The relation between the velocity derived from the X-ray temperature and from the 
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FIG. 1: Velocity statistics for all clusters larger than a minimum mass Mmin, evaluated for two values of Mmin differing by 
20%. (a) The mean pairwise streaming velocity and (b) the perpendicular velocity correlation function. 



electron temperature is just 

?^truo = {a/Tx + 6)Wobs, (1) 

where Utruc is the velocity inferred from the electron temperature and Uobs is the velocity inferred from the X-ray 
temperature; Sunyaev-Zeldovich cluster velocity measurements will be correct when using the electron temperature. 
For a cluster of temperature Tx = 3 keV, the first term is around 10% of the second term, and the relative contribution 
decreases further for more massive clusters; we thus neglect the first term, leaving the true velocity proportional to 
the observed velocity, 

Wobs = Pvtrnc (2) 

with (3 = 1/h. For the given values oi h ~ 0.69 and 0.53, the cluster velocity derived using a measured X-ray 
temperature will be 1.4 and 1.9 times larger than the true velocity derived from the gas temperature. To the extent 
that we will not know perfectly the Tx-Tg relation, our inferred cluster velocities will be dominated by our fractional 
misestimate of (3: a 10% overestimate of (3 gives about a 10% overestimate of the cluster velocity. Note that for 
a sample of clusters, b (and thus /3) will likely be easier to infer than a since it represents the slope of the T^-Tx 
relation, rather than an extrapolation of this relation to Tx = 0. 

Point sources can be modeled in the same way; Ref. [131 discusses the fact that systematic errors in velocity from 
point sources can be significant even though their contribution to statistical error may be small. We again assume 
an observed velocity proportional to the true velocity; we also add a constant offset, Vtmc = Pvohs — VoS- Numerical 
simulations with large numbers of clusters suggest the value /3 — 2, perhaps slightly larger than that expected from 
X-ray temperatures. 

To model these effects, we consider a hypothetical cluster velocity survey with velocities measured for all clusters of 
mass M > 2 X 10^'^h~^ Mq with a statistical error of dy = 300 km/s, plus cosmic variance and Poisson noise for a given 
survey area. This is a rather optimistic scenario in terms of the measurement errors. We are interested in quantifying 
the degradation in dark energy constraints due to systematic errors, compared to the case of only statistical errors. 
The most optimistic statistical error gives an estimate with maximum relative systematic error impact; hence, the 
results here represents the worst case scenario. We also assume a rather high value for erg — 0.9. For a lower value 
as ~ 0.75, the total number of clusters above a given mass threshold will be reduced by 40%, which in turn will 
increase the statistical errors accordingly. 

Fig. [2] shows the effect of systematic errors on the pairwise mean streaming velocity and the perpendicular velocity 
correlation function, as determined from the lightcone output of the VIRGO simulation [2^. The dashed line sur- 
rounded by the shaded region shows the actual value of the statistics as drawn from the simulation, with inferred 
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statistical errors assuming a velocity error of for a 4000 deg^ sky area; see Ref. [Tj] for details. The higher offset solid 
lines show the same quantities except with the individual cluster velocities biased using /3 = 1.7 and /3 = 2, while the 
dot-dash line shows a constant velocity offset corresponding to Uofi = 30 km/s. Note that a constant offset has no 
effect on the mean pairwise streaming velocity, and a relatively small effect on the correlation function compared to 
the shift due to (3. Constant velocity offsets would also be evident in the velocity distribution function, since the entire 
cluster peculiar veloctiy sample should have zero mean [l3l |. For the mean pairwise velocity, a bias corresponding to 
/3 = 1.7 shifts the velocity statistic by about la statistical error, while the effect is substantially larger in the velocity 
correlation function. 

In the case of X-ray temperature, we already have reasonable estimates of the difference in X-ray and electron 
temperatures, from both analytic and numeric calculations; the actual bias in velocity statistics will be due only to 
our error in understanding this relation, which should be much smaller than the size of the effect displayed in Fig. [51 
The extent to which we can characterize and understand the effect of the point source population is currently under 
investigation and requires a better observational characterization of the relevant sources and their correlation with 
galaxy clusters. Ultimately the point sources can be spatially resolved by observations at submillimeter wavelengths, 
but doing this over a survey region of several hundred square degrees is likely impractical in the foreseeable future. 

The numbers presented here are a worst-case scenario, should we have a gross misunderstanding of cluster physics, 
or completely fail to recognize a substantial source of systematic error in cluster velocity estimates. We have simply 
assumed that we do not account for systematic offsets in X-ray temperature compared to electron temperature, or 
systematic errors in cluster velocity estimates due to point source contamination. We already have detailed estimates 
of the former, based on simulations, and the latter is under active study [H, |2^, [s^, HH . We can also hope to measure 
these systematic effects directly from the cluster velocity data; such "self-calibration" will be considered below. 



IV. BIAS IN DARK ENERGY PARAMETERS 



In order to study the bias induced in dark energy parameters from systematic errors in the velocity statistics, we 
consider a fiducial cosmology described by the set of cosmological parameters p on which the velocity field depends: 
the normalization of the matter power spectrum erg, the power law index of the primordial power spectrum ns, and 
the Hubble parameter ft., plus the dark energy parameters namely, JIa and two parameters wq and Wa describing 
the redshift evolution of its equation of state w(a) = wq -\- {1 — a)wa- We assume Gaussian priors with variances of 
AcTg = 0.09, Arts — 0.015 [Hj and Aft — 0.08 [33|. We then perform a simple Fisher matrix analysis to find the bias 
on these parameters from measurements of our two velocity statistics with small systematic errors. 

The Fisher information matrix for each of the two statistics is [l2| 

^.. = Ef^[^ffe-)r^f^ (3) 

where stands for either Vij{r, z) or {viVj){r, z), C"^(ij) is the total covariance matrix in each bin. A detail description 
of the statistics and its covariance calculation are given in Ref. [l3|- Assuming the systematic offsets in the velocity 
statistics are small so that the gaussian assumption is valid, the bias in parameter p can be written as [s^l 




where (/)sys — Scj) the difference between the biased and the true value. Note that the assumption of small offsets may 
not be valid in our case; nevertheless it gives us an estimate of the magnitude of the bias. 

Tables U & |ll] show the dark energy parameter biases for each of the two velocity statistics, and for each of two 
survey areas. Assuming a measurement error normally distributed with — 300 km/s, both Vij and (viVj) put tight 
constraints on dark the energy density and relatively weak constraints on its equation of state in the absence of any 
systematic bias. The systematic bias for wo and Wa is only marginally greater than the no-bias statistical error for 
(3 = 2, for both the statistics and the survey areas considered; the only exception is for the case of Vij for a 4000 deg^ 
survey, where the systematic bias on Wa is 3 times greater than the statistical error. For 17 a the bias is substantial 
for all survey areas and both the statistics; Vij generally gives a smaller bias on fl\ than {viVj). While f^A is strongly 
constrained by other measurements, one virtue of a velocity survey is a completely independent constraint on f2A- 
Introducing prior cosmological constraints consistent with projections for the Planck satellite, velocity statistics will 
also provide competitive constraints on wq and Wa [l3j. Hence it is important to determine whether self-calibration 
of unknown systematic errors will help reduce the bias in determining these parameters. 

We note in passing that non-trivial covariances exist between the various parameters; these will be considered in 
greater detail elsewhere using Monte Carlo explorations of the parameter space. 
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Survey Area 


Parameter (p) la 


error (Ap) bias (Sp) 


5p/Ap bias (Sp) 


Sp/Ap bias (Sp) 


Sp/Ap 


deg2 






P = 1.4 


P = 1.4 


P = 1.7 


P = 1.7 


P = 2.0 


P = 2.0 




nA [0.7] 


0.03 


0.04 


1.46 


0.075 


2.56 


0.11 


3.7 


4000 


wo [-1] 


0.35 


0.15 


0.43 


0.35 


0.75 


0.37 


1.1 




Wa [0] 


0.56 


0.65 


1.15 


1.13 


2.0 


1.6 


2.9 




[0.7] 


0.034 


0.07 


2.0 


0.12 


3.5 


0.17 


5.1 


400 


wo [—1] 


0.84 


0.07 


0.08 


0.12 


0.14 


0.17 


0.2 




Wa [0] 


1.5 


0.7 


0.5 


1.24 


0.8 


1.8 


1.2 



TABLE L The statistical errors Ap in dark energy parameters Qa, Wq, and Wa, and the bias Sp in these parameters due to 
systematic error in cluster velocity estimates, using the mean pairwise streaming velocity Vij{r). The fiducial cosmological 
model has Us = 1, erg = 0.9, h = 0.7, 57a = 0.7, wo = —1, and Wa = 0, with prior normal errors of Aus = 0.015, Acrg = 0.09 and 
Ah = 0.08 and a spatially flat universe assumed. No priors on dark energy parameters are included. Cluster velocity normal 
errors of = 300 km/s are assumed. 



Survey Area Parameter (p) la 


error (Ap) bias (Sp) 


Sp/Ap bias (Sp) 


Sp/Ap bias {Sp) 


Sp/Ap 


deg2 






P = 1.4 


P = 1.4 


P = 1.7 


P = 1.7 


P = 2.0 


P = 2.0 




S^A [0.7] 


0.038 


0.22 


5.7 


0.43 


11.5 


0.68 


17.8 


4000 


wo [-1] 


0.41 


0.09 


0.22 


0.17 


0.43 


0.28 


0.69 




Wa [0] 


0.71 


0.31 


0.43 


0.6 


0.85 


0.96 


1.35 




^A [0.7] 


0.08 


0.2 


2.7 


0.4 


5.3 


0.7 


8.4 


400 


wo [-1] 


0.96 


0.11 


0.12 


0.22 


0.23 


0.35 


0.36 




Wa [0] 


1.9 


0.5 


0.27 


1.01 


0.54 


1.6 


0.86 



TABLE II: The same as in Table U but for the velocity correlation function. 




FIG. 2: The efi'ect of systematic errors for the mean pairwise streaming velocity (left) and the velocity correlation function 
(right). The shaded region with the dashed line shows each statistic obtained from the Virgo dark-matter simulation. The 
effect of unmodeled systematic bias between estimated and actual gas temperature are shown with the red solid line (/3 = 1.7) 
and the blue dashed line (/3 = 2). Also shown as dot-dashed black lines are the effects of a constant velocity offset VoB = 30 
km/sec. Note that Vij has the advantage of being insensitive to velocity offsets. 
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Survey Area Parameter (p) la error (Ap) self-cal(Asp) Asp/Ap self-cal(Aap) Asp/Ap self-cal(Asp) Asp/Ap 



deg2 






13= 1.4 


13= 1.4 


13 = 1.7 


13 = 1.7 


P = 2.0 


13 = 2.0 




f^A [0.7] 


0.029 


0.03 


1.03 


0.031 


1.07 


0.032 


1.10 


4000 


wo [-1] 


0.35 


0.37 


1.06 


0.38 


1.08 


0.41 


1.175 




Wa [0] 


0.56 


0.71 


1.27 


0.75 


1.35 


0.78 


1.4 




[0.7] 


0.034 


0.036 


1.06 


0.036 


1.06 


0.037 


1.07 


400 


wo [—1] 


0.84 


0.87 


1.04 


0.89 


1.06 


0.93 


1.11 




Wa [0] 


1.5 


1.7 


1.13 


1.9 


1.27 


2.2 


1.46 



TABLE III: Constraints with the self-cahbration of the systematic parameters /3o and 7 for Vij. A prior of ±0.5 is assumed on 
the systematic parameters namely /3o & 7- Note that AsP & Ap denotes the 1 — a statistical error on dark energy parameters 
when nuisance parameters are included and not included respectively 



Survey Area Parameter (p) la error (Ap) self-cal(Asp) Asp/Ap self-cal(Asp) Asp/Ap self-cal(Asp) Asp/Ap 



deg^ 






P = 1.4 


P= 1.4 


P = 1.7 


P = 1.7 


P = 2.0 


P = 2.0 




^A [0.7] 


0.038 


0.042 


1.09 


0.044 


1.15 


0.047 


1.23 


4000 


wo [-1] 


0.41 


0.65 


1.6 


0.8 


1.97 


1.0 


2.4 




Wa [0] 


0.71 


2.64 


3.7 


3.1 


4.4 


3.6 


5.2 




^A [0.7] 


0.08 


0.088 


1.09 


0.095 


1.18 


0.11 


1.31 


400 


wo [-1] 


0.96 


1.6 


1.69 


1.97 


2.0 


2.42 


2.52 




Wa [0] 


1.9 


5.3 


2.8 


6.7 


3.6 


8.3 


4.4 



TABLE IV: Same as Table |III] but for {viVj). 



V. SELF-CALIBRATION OF SYSTEMATIC OBSERVABLES 

One potential method for dealing with systematic errors is to adopt some reasonable parameterized model for the 
errors, then solve for these systematic error parameters along with the cosmological parameters of interest, given 
the data in hand. In this section we study the self-calibration of the variable that models systematic velocity errors 
due to X-ray temperature offset or imperfect point-source subtraction. To this end, we allow (3 to co-vary with the 
cosmological parameters. To allow for redshift evolution of /?, we write (3 = /3o/(l -f z)"' with [3o being the value 
of /3 at z = 0. We choose a fiducial value of 7 = 0, while fiducial values of (3q — 1.4, 1.7, and 2 correspond to 
the three values of (3 considered in Sec. IIVI We also assume a mild normal distribution prior on these parameters 
with a variance of 50% . We envisage such moderate priors can be obtained using numerical simulation studies and 
follow-up observations. We perform a fisher matrix analysis using Eq. ([3]) with the cosmological parameters plus the 
two systematic parameters /3o and 7, then marginalize over /3o and 7 to get the dark energy parameter constraints. 
Cosmological parameters will necessarily have their constrains weakened, but if the model for the systematic errors 
is an accurate representation of the actual systematic errors, the bias in cosmological parameters will be reduced. 

Results are shown in Table IIIII & IIVI for the two statistics and (viVj) and two survey areas. For Vij., the 
degradation in the constraint on JIa varies from 3% to 10% for f3 = 1.4 to 2.0. For wq, the degradation varies from 
6% to 18% and for Wa it is 27% to 40% for a 4000 deg^ survey area. For 400 deg^ the relative degradation of the 
constraints is smaller since the statistical error is larger than for greater sky area. For (uiVj) the degradation is larger 
than for Vij since it varies as . Table |V] gives the constraints on the parameters /Sq and 7 that are used to describe 
the systematic velocity errors: Vij gives 20% and 30% constraints on /3q and 20% and 40% constraint on 7 for the two 
survey areas 4000 and 400 deg^ respectively. 

VI. SUMMARY AND CONTEXT 

Galaxy cluster velocity surveys, with velocities determined from kinematic Sunyaev-Zeldovich measurements, have 
the potential to provide constraints on dark energy parameters competitive with other methods. Their challenge is 
attaining sufficient sensitivity over a large enough sky region; the kinematic Sunyaev-Zeldovich signal is typically only 
a tenth of the thermal SZ signal for the same galaxy cluster. But while cluster number counts based on the thermal SZ 
signal probe dark energy with a signal that is easier to detect, cluster velocities have the tangible advantage that they 
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Survey Area Vij {viVj) 

deg2 /? = 1.4 = 1.7 /3 = 2.0 = 1.4 /3 = 1.7 /3 = 2.0 

4000 0.24 [0.25] 0.21 [0.23] 0.18 [0.22] 0.09 [0.2] 0.085 [0.16] 0.08 [0.13] 

400 0.32 [0.42] 0.3 [0.4] 0.27 [0.39] 0.16 [0.38] 0.17 [0.33] 0.17 [0.29] 

TABLE V: Constraints on the parameters /3o and 7 used to model the systematic offset. Constraints are shown as A/3o [A7]. 

are nearly insensitive to the uncertainties in the connection between the underlying cluster mass and the observed 
cluster SZ signal that will substantially impact cluster number counts. We have demonstrated this fact explicitly by 
computing the cluster mean streaming velocity and velocity correlation function from the halo model, finding that 
both statistics are nearly independent of the cluster mass threshold chosen. Galaxy clusters serve as tracers of the 
underlying large-scale velocity field, independent of their masses. 

Cluster velocities are not without their own systematic errors, however; in particular, the velocities inferred from 
kSZ measurements may have systematic errors arising from a mis-estimated relation between the cluster X-ray tem- 
perature and its gas temperature, or from an incorrect correlation between the thermal SZ signal and the cluster 
gas temperature. Either might be used to extract the cluster velocity from three-band microwave measurements. 
Further complications arise from infrared point sources, which have variable spectral indices and can be substantially 
correlated with galaxy cluster positions. We have considered a simple toy model which assumes a linear relationship 
between the measured and actual cluster velocity. If this relation is not well understood, we have quantified the bias it 
can induce in inferred values of cosmological parameters. This bias can be significant compared to the corresponding 
statistical errors, although reasonable levels of understanding of these correlations will likely reduce the systematic 
errors to below the level of statistical errors. 

This work has considered the cluster peculiar velocity to be the basic observable. Velocity statistics are relatively 
straightforward to compute in particular cosmological models [ij, [13, [iM [3 1 t)ut the kinematic SZ effect actually 
measures a line-of-sight integral of the gas density times the peculiar velocity, proportional to the cluster gas mo- 
mentum. The systematic errors associated with inferring the cluster velocity, particularly the discrepancy between 
X-ray temperature and gas temperature, can be avoided completely by measuring the cluster gas momentum via the 
kinematic SZ distortion, then comparing these to theoretical distributions of cluster momenta instead of velocities. 
The difficulty then becomes accurate prediction of cluster momentum statistics, which must be done from simulations. 
While total cluster momentum is relatively straightforward to extract from large-volume cosmological dark matter 
simulations, the cluster gas momentum depends on the cluster gas fraction. This fraction is expected on general 
grounds to be relatively constant for large clusters, since it is difhcult for gas to escape from their deep gravitational 
wells, but their star formation history might vary significantly and depend in detail on the physics of star formation 
and energy feedback from stars and active galaxies. This question is addressed for lower-mass galaxy groups in [sH, 
but realistic estimates for galaxy clusters await larger- volume cosmological hydrodynamic simulations. Uncertainty 
in cluster baryon fraction could be a significant systematic error when using cluster momentum statistics to constrain 
cosmology. 

To date, a galaxy cluster peculiar velocity has never been measured directly: no Sunyaev-Zeldovich measurements 
have had the combination of arcminute angular resolution and micro-Kelvin sensitivity needed to make a statistically 
significant detection of the kinematic SZ signal. (Velocity upp er Hmits have been obtained for a few clusters [SGllSTj.l 
This situation will change in the near future; both ACT [i|,[3|] and SPT likely will have sufficient angular resolution, 
sensitivity, and frequency coverge, and are currently taking data. The first direct measurement of the peculiar velocities 
of objects at cosmological distances will be a signature achievement, and when extended to large surveys will give 
an independent route to constraining dark energy properties [Tsj and properties of gravitation (S. Bhattacharya & 
A. Kosowsky, in preparation (2008)). Capitalizing on this opportunity requires sufficient control over systematic 
errors, as with all dark energy probes. As we have shown here, cluster velocities have a large advantage over cluster 
number counts because the selection function of the clusters has little impact on the cluster statistics we want to 
measure; cluster number counts, in contrast, are equally sensitive to the selection function as to cosmology. Systematic 
errors in velocity estimates for individual clusters can potentially cause significant bias in the inferred values of dark 
energy parameters, but based on the calculations presented here we anticipate that careful observations and modeling 
of galaxy clusters can limit these systematics to levels below the statistical errors in cosmological constraints. 
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